Hadroproduction of the xi an d X2 States of Charmonium 
in 800 GeV/c Proton-Silicon Interactions 
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The cross sections for the hadroproduction of the Xi an d Xi states of charmonium in proton-silicon collisions at 
y/s = 38.8 GeV have been measured in Fermilab fixed target Experiment 771. The \ states were observed via their 
radiative decay to J/ipj, where the photon converted to e + e~ in the material of the spectrometer. The measured 
values for the xi an d X2 cross sections for xf > are 263 ± 69(stat) ± 32(syst) and 498 ± 143(stat) ± 67(syst) nb 
per nucleon respectively. The resulting cr(xi)/&{x2) ratio of 0.53 ± 0.20(stat) ± 0.07(syst), although somewhat larger 
than most theoretical expectations, can be accomodated by the latest theoretical estimates. 
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PACS numbers: 13.85.Ni, 13.85. Qk, 14.40. Gx, 25.40.Ve 

Charmonium hadroproduction has provided interest- 
ing challenges to the understanding of QCD. Early at- 
tempts to describe the formation of a cc bound state, 
according to the color evaporation jlj or color singlet Q 
models, did not provide a satisfactory description of the 
available data. More recently, the Non-Relativistic QCD 
Factorization Approach j^] , which incorporates in a more 
rigorous fashion some of the features of the previous mod- 
els, has provided a more successful description of the pro- 
cess. For any model, a rather crucial test has been the 
prediction of the relative rate of production for differ- 
ent charmonium states. In particular, when dealing with 
proton-induced processes, the absence of quark annihila- 
tion diagrams and the suppression in Leading Order of xi 
production gluonic diagrams implies rather small values, 
typically less than 10%, for the ratio of xi 1° X2 produc- 
tion The measurement presented here represents a 
significant contribution to the available data, since it is 
the first observation of cleanly resolved %i and xi states 



in a proton- induced fixed target experiment. 

The FNAL E771 experiment utilized a large- 
acceptance spectrometer Q to measure several pro- 
cesses containing muons in the final state. Protons of 
800 GeV/c momentum were transported by the Fermilab 
Proton West beam line to the High Intensity Laboratory, 
where they hit a 24 mm thick silicon target. Operating 
at a beam intensity of ps 3.6 x 10 7 protons per spill- 
second, the experiment accumulated a total of 6.4 x 10 11 
p-Si interactions. The incoming proton beam trajectory 
and flux were measured by a six plane silicon detector 
station. The 0.26 radiation length target was composed 
of twelve 2 mm silicon foils separated by 4 mm. The tar- 
get was followed by a microvertex detector consisting of 
fourteen 300 /im thick silicon planes that, while not used 
in the analysis presented here, contributed an additional 
0.045 Xo to the target region radiation length. 

The spectrometer's tracking system consisted of seven 
multi-wire proportional chambers and three drift cham- 
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bers upstream plus three drift chambers and six combi- 
nation drift /pad/strip chambers downstream of a dipole 
analysis magnet which provided an 821 MeV/c pt kick in 
the horizontal plane. Downstream of the wire chamber 
system, an electromagnetic calorimeter consisting of an 
active converter and 396 scintillating glass and lead glass 
blocks was used for electron/positron identification. The 
final clement of the spectrometer, a set of three planes 
of resistive plate counters (RPC's) segmented into 512 
readout pads and sandwiched between layers of steel and 
concrete absorbers, provided muon identification. The 
material in the absorber walls represented an energy loss 
of 10 GeV in the central region and 6 GeV in the outer 
region of the detector for the incident muons. 

A dimuon trigger || selected events with a J/ip in the 
final state via the decay J/ip — » /i + /J~. A trigger muon 
was defined as the triple coincidence of the OR of 2 x 2 
pads in the first RPC plane and the OR of 6 x 6 pads in 
the second and third RPC planes in projective arrange- 
ments. A dimuon trigger was defined as two such triple 
coincidences. The trigger reduced the 1.9 MHz interac- 
tion rate by a factor of ~ 10 4 , selecting approximately 
1.3 x 10 8 dimuon events to be written to tape. 

The seed for muon track reconstruction was provided 
by the RPC triple coincidences. The roads formed by 
the pads involved in the coincidences were projected into 
the rear chamber set, identifying a region in which to 
search for candidate muon tracks. Muon tracks recon- 
structed in the rear chamber set were then matched with 
tracks found in the front chamber set by requiring a good 
front-rear linking \ 2 - Muon pairs were required to come 
from a common vertex by applying a cut to their dis- 
tance of closest approach. About fifty thousand dimuon 
events survived the reconstruction process, quality cuts, 
and vertex cuts. 
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FIG. 1. Opposite sign dimuon mass spectrum. The solid 
line is a fit to the signal plus background and the dashed line 
is a fit to the background only. Inset shows the T signal. 

Figure 1 shows the resulting dimuon mass spectrum 
containing peaks corresponding to the J/ip, ip(2S), and 



T (inset). Superimposed on the dimuon mass spectrum is 
a fit to the data obtained with the sum of two Gaussians 
for the J/ip peak, a single Gaussian for the ip(2S), and 
the form — §— exp{— 6m„„) for the continuum background. 

The two Gaussians fit to the J/ip peak is a good approxi- 
mation (as confirmed by Monte Carlo) to a non-constant 
mass resolution, caused by the confusion associated with 
increases in hit density near the beam region. The num- 
ber of J/-0's and ip(2S)'s after background subtraction 
was 11,660 ± 139 and 218 ± 24, respectively §. 

Events in a window of ± 100 MeV/c 2 around the J/ip 
mass were refit varying the muon momenta within mea- 
surement errors, with the constraint that the invariant 
mass of the pair be equal to J/ip mass. The result- 
ing dimuon event sample was then inspected to search 
for e + e~ pairs that might be the result of conversions 
of photons from x ~ * J/ip 7 decays. Dimuon events 
which contained pairs of tracks matching the topology 
of a 7 — > e + e~ conversion in the target region - collincar 
before the magnet in both bend and non-bend projec- 
tions, collinear in the non-bend plane and coplanar in 
the bend plane after the magnet - were then designated 
as x decay candidates. 

All electron/positron pair candidates were required 
to satisfy additional conditions. At least one of the two 
track candidates was required to be associated with an 
energy deposition in the calorimeter consistent with an 
electro-magnetic shower. In addition, the total trans- 
verse momentum of the e + e~ pair in the rear of the mag- 
net was required to be zero (within the resolution of the 
spectrometer) relative to the common e + e~ trajectory in 
front of the magnet. To quantify how well a pair fitted 
the 7 — * e + e~ hypothesis, a x 2 wa s formed, 

2 ( a xl - a x2) 2 (Oyl - a y 2) 2 (fiyl ~ b y2 ) 2 , . 

er 2 + er 2 cr 2 +cr 2 er 2 + cr 2 [ 1 

u ax\^ u ax2 " ayl T u ayl ° byl ~ ° by2 

where a x i and a x2 are the electron and positron track 
intercepts at the magnet in the bend plane, a v \ and a y2 
the track intercepts in the non-bend plane, b y \ and b y2 
the track slopes in the non-bend plane, and the er's are 
the measurement errors on these quantities. The elec- 
tron/positron candidate with the smallest x 2 m a given 
event was designated as a photon conversion candidate. 
Additional cuts requiring a good x* \ the transverse mo- 
mentum of the parent photon to be between 250 and 700 
MeV/c, and the invariant mass squared of the e + e~ pair 
to be less than 3000 (MeV/c 2 ) 2 were applied to maxi- 
mize signal to background in the final sample of events 
containing both a J/ip and a photon conversion. 

The J/ip e + e~~ invariant mass shown in Fig. 2 was 
calculated using the electron and positron momenta ob- 
tained from the tracking system. Clear xi an< i X2 sig- 
nals can be seen. The background to the xi an d X2 was 
well-described by uncorrelated e + e _ and J/ip combina- 
tions: the solid line of Fig. 2 was obtained by fitting 
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two Gaussians plus a polynomial background. The poly- 
nomial background was obtained by fitting the mass dis- 
tributions of J/ip's and e + e~'s extracted from different 
events. The numbers of xi and X2 obtained from the fit 
are 33 ± 9 and 33 ± 10, respectively. The fitted width is 
5.2 ± 2.0 MeV/c 2 for both the xi and X2 peaks. 
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FIG. 2. J/tpe + e~ invariant mass. The solid line shows a 
fit to a polynomial background plus two Gaussians of equal 
width. 



an e+e - identity. This sample was then subjected to 
the same constraints as those applied in the x analysis. 
The resulting overall acceptance times efficiency (inclu- 
sive of conversion probability) for photons from x decay 
was determined to be (8.25 ± 0.4) x 10~ 3 . 

Using the 7 — > e + e~ acceptance and efficiency, the 
measured branching ratios for xi and X2 into J/ip 7 [jl0| , 
the measured J/ip pN forward cross section at 
V / s=38.8 GeV and the number of observed xi> X2 and 
J/ip — ► the absolute xi and X2 cross sections for 

xp > were calculated to be er(xi) = 263 ± 69(stat) ± 
32 (syst) nb/nucleon and cr(x 2 ) = 498 ± 143(stat) ± 67 
(syst) nb/nucleon, respectively. The main contributions 
to the systematic errors came from the error on the J/ip 
cross section (9%), the uncertainty in the knowledge of 
the cut efficiencies (5%) and the errors on the branching 
ratios for xi (6%) and X2 (8%) . 

Using the production cross sections for xi and X2, 
the ratio of the xi to X2 production cross sections was 
determined to be ct(xi)/c(x2) = 0.53 ± 0.20(stat) ± 
0.07(syst). Combining this result with the two previous 
measurements of x production by a proton beam [pi, 
we have computed the world average (shown in Fig. 3) 
to be cr(xi)/ (J (X2) = 0.31 ± 0.14. This figure is consistent 
with the latest NRQCD estimates of w0.3 where xi 
production was boosted by the inclusion of higher order 
terms in the velocity expansion G3]. 



To determine the total cross section for xi and X2 
production, the overall acceptance times efficiency for 
photon conversion and for electron/positron acceptance 
and recontruction efficiency had to be determined. To 
accomplish this, a Monte Carlo sample of % — ¥ J/ip 7; 
J/ip — > decays was generated using Pythia ||]. The 

photon and the muons were then propagated through a 
GEANT H simulation of the E771 detector, including 7 
conversion, scattering, bremsstrahlung and dE/dx. Hits 
from the Monte Carlo tracks obtained by this prescrip- 
tion were then inserted into actual dimuon trigger events 
to simulate realistically backgrounds and losses in pattern 
recognition due to confusion from noise hits and other 
tracks. Measured detector efficiencies were also applied 
to the inserted hits. These hybrid Monte Carlo and data 
events were analyzed in a manner identical to the data in 
order to determine acceptances and tracking efficiencies. 

Rather than attempting to simulate the electromag- 
netic calorimeter response to e ± in detail in a Monte 
Carlo, the efficiency of matching an electron or positron 
candidate to a shower in the calorimeter was determined 
using a large sample of electron/positron pairs from pho- 
ton conversions in minimum bias events. A sample of 
e + e~ pairs with kinematics similar to those of the x 
e + e~ pairs was collected using very tight cuts to ensure 
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FIG. 3. o"(xi)/°"(X2) in proton-nucleon interactions. Re- 
sults from this experiment (E771) are shown, together with 
the two existing results and the world average with la errors. 

Finally, the energy dependence of the combined xi 
and X2 production near threshold was compared to the 
corresponding quantity for J/ip production. In Ref. |TJ] 
data on J / %p production from seventeen different pN ex- 
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periments over a large range of center-of-mass energy, 
y/s s; 8 to 52 GeV, were fit as a function of y/s. The 
J/ip production data near threshold was well represented 
by the function <j{\fs)j/^ = cr (l — Mjj^/^fs)®, with 
cr = 1.0 ± 0.1 ^b/nuclcon and (3 = 11.8 ± 0.5. To check 
whether \ production has similar dynamics as J/ip pro- 
duction, the sum of the x cross sections have been fit 
to a similar parameterization with (3 fixed to the J /ip 
value and M x replacing Mj/^,. The result of the fit, 
shown in Fig. 4, shows the similarity of the J/ip thresh- 
old production parametrization to the threshold behav- 
ior of the combined x state cross sections. The fit yields 
(j o = 2.3 ± 0.4/ib/nucleon for the asymptotic <j(x) cross 
section. 
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FIG. 4. a(x) vs. \fs in proton-nucleon interactions. Su- 
perimposed on the data is a fit to a threshold production 
parameterization described in the text; here <r(x) is the sum 
of the xi and X2 cross sections where the E673 data point has 
been obtained from the published xi cross section and ratio 
of X2 to xi 
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